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ABSTRACT 

Ultraluminous X-ray sources (ULXs) with X-ray luminosities larger than the Eddington lu- 
minosity of stellar mass objects may be powered by intermediate mass black holes (IBHs) 
of masses M. ~ 10 3 M Q . If IBHs form in young dense stellar clusters, they can be fed by 
Roche lobe overflow from a tidally captured massive (M* > 10Af Q ) stellar companion. Af- 
ter the donor leaves the main sequence it forms a compact remnant, which spirals in due to 
gravitational wave (GW) emission. We show that space based detectors such as the Laser 
Interferometer Space Antenna are likely to detect several of these sources. GW sources stem- 
ming from this scenario have small eccentricities which give distinct GW signals. Detection of 
such a GW signal will unambiguously prove the existence of IBHs, and support the hypothesis 
that some ULXs are powered by IBHs with captured companions. 

Key words: gravitational waves - black hole physics - stellar dynamics - galaxies: star clusters 
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1 INTRODUCTION 

The counterpart of ultraluminous X-ray sources (ULXs), which 
have X-ray luminosities Lx larger than the Eddington luminos- 
ity of a stellar mass object of mass M, Lx > Le = 1.3 x 
10 39 erg s _1 M/WMq, are not known. Most likely there is no uni- 
versal engine for ULXs . Some may be powered by anisotropic ra- 
diation (e.g. King et al j 200 ll : Rappaport, Podsiadlowski & Pfahl 
2005) or super Eddington luminosity (Begelma rl 20021) from stel- 
lar mass black holes. 

A third possibility is offered by the hypothesised existence of 
intermediate mass black holes (IBHs; see Miller & Comer 2004 for 
a review) of masses 10 3 M Q < M. < W 5 M Q , which could radi- 
ate isotropically at a sub-Eddington rate to account for the observed 
luminosities. There has been some observational evidence that at 
least some ULXs are powered by IBHs. In some cases the mass- 
tempe rature relation (Miller et al j 20031 Miller, Fabia n & M iller 
I 2004[) or q uasi periodic oscillations (Fiorito & Titarchuk 200^ Liu 
et al 2005) indicate a high accreting mass. 

Numerical N-body and Monte Carlo simulations indicate that 
the stellar collision rate in young dense stellar clusters during core 
collapse can become very large, giving rise to a hierarchical merger. 
In that case an object of thousands of solar masses may form (Porte- 
gies Zwart et al l 199911 2004 Giirkan, Freitag & Rasi d 20041: Fre - 
itag, Giirkan & Rasic 2005a; Freitag, Rasio, & Baumgardi 2005b). 
The fate of such a massive star is unclear, but it might lead to the 
formation of an IBH. This scenario is supported by the f act that 
ULXs correlate positively with star formation (Swartz et al l 20041 ; 
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Liu, Bregman & Irwi nl 20051) and t hat som e ULXs are associated 
with stellar clusters (e.g. Zezas et al l 20021) . 

The young stars in the host cluster of the IBH have strong 
winds which blow out the gas from the cluster, and there is insuf- 
ficient free gas available to power a ULX. However, the IBH can 
acquire a companion star by dynamical capture or tidal capture. 
Here we discuss the latter possibility. Hopman, Portegies Zwart & 
Alexander (2004) showed that tidal capture of a main sequence star 
of mass Mi, and radius i?* can lead to circularization close to the 
tidal radius 

n - «.(£)"', 

which is the distance from the IBH where the tidal forces equal the 
forces which keep the star bound. As the star evolves it starts to fill 
its Roche lobe causing it to lose mass to the IBH. This leads to high 
X-ray luminosities, provided that the donor is sufficiently massive 
(M^ > I QMp, H opman et al. l 20041: Portegies Zwart, Dewi, & 
Maccarone 2004f L i 20041) . Interestingly, for some ULXs an opti- 
cal counterpart has been identified, indicating that these ULXs are 
binary systems (Liu, Bregman, & Seitze J 2004 Kuntz et a l l 20051) . 

In this scenario a ULX may turn on for at most the life- 
time t* ~ 10 Myr of the captured star. In addition to strong X- 
ray emission the star emits gravitational waves (GWs). Portegies 
Zwart (2004) discussed the possibility of observing X-rays and 
GWs simultaneously. Only if the captured star is sufficiently light 
(Mi, < 2M©), its tidal radius is small enough for it to emit GWs 
in the LISA band during the mass transfer phase (Portegies Zwart 
I 2004 . However, in order to account for the high ULX luminosi- 
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ties the companion star should be considerably more massive than 
2 M Q (Hopman et al l 2004 Portegies Zwart et al l 2004 L i 20041) . 
In young and dense star clusters mass segregation causes most mas- 
sive stars to accumulate in the cluster centre, and the combination 
of a high mass and large size make these stars excellent candidates 
for tidal capture. When after time t* the massive donor explodes it 
turns into a compact remnant (CR); a neutron star (NS) or a stellar 
mass black hole (SBH). From that moment the ULX, deprived of 
its source of gas, turns off. 

Here we focus on the subsequent -post supernova- evolution 
of the (IBH, CR) binary. In-spiral of the CR into an IBH due to 
GW emission results in a strong signal in frequencies measurable 
by space detectors such as LISA, providing a wealth of information 
about the system (see Mille J2OO2I for a discussion). We show that 
LISA may be able to observe such (IBH, CR) binaries. 



2 STELLAR CAPTURE AND BINARY EVOLUTION 

When stars orbit IBHs on wide, but highly eccentric orbits, with 
periapse close to the tidal radius of the IBH, r p > r±, orbital en- 
ergy is invested in tidal d istortions of the star, causing it to spiral 
in. Hopman et al. 1(20041) showed that if the IBH is less massive 
than Af max ~ 1O 5 M , the star may be able to cool down by ra- 
diating the excess energy efficiently, and survive the strong tidal 
forces. Eventually it circularises near the IBH. We first estimate the 
fraction /merge of isolated binaries which merge as a result of GW 
emission within the age of the Universe, while accounting for an- 
gular momentum conservation during the mass transfer phase and 
an isotropic velocity kick. We then discuss the consequences of the 
gravitational interaction of cluster stars with the (IBH, CR) binary. 



2.1 Evolution of an isolated IBH-star binary 

After tidal circularization near the IBH, the star may fill its Roche 
lobe, possibly after a period of stellar evolution and expansion. 
Roche lobe overflow from a A/* > 1OM donor can then give 
rise to high lumino sities f or a time limited by t he life-t ime of the 
star (Hopman et al l 20041 ; Portegies Zwart et al l 2004 l A 2004 . 

The main-sequence star fills its Roche lobe at a distance ~ 2r t 
before the terminal age main sequence. We assume that the entire 
hydrogen envelope is transferred to the IBH, while the binary con- 
serves angular momentum. The mass of the stellar core of the donor 
star is given by M c = 0.08 Mq(M+/Mq) 1a (Iben, Tutukov 
& Yungelson|j993), while the mass of the hydrogen envelope is 
Mh = — M c . Angular momentum conservation implies that 
the binary separation increases during the mass transfer phase. 

The remaining helium star subsequently explodes in a su- 
per nova after a time t*, and forms a CR. In this event a star of 
1OM < M t < 2OM forms a NS of Chandrasekhar mass, 
Mch = 1.4M . Stars more massive stars than 20 M collapse to 
SBHs; we assume that the mass distribution is that found by Fyer 
&Kalogera lf200ll) . 

The semi-major axes and eccentricities of the CR are deter- 
mined by the post mass-transfer orbital elements while accounting 
for the mass lost in the super nova and for the velocity kick, im- 
parted to the CR at the time of the explosion. The kick velocity 
is taken in a random direction with magnitude from the distribu- 
tion of pulsars velocities, which is well fitted by a double Gaussian 
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where w = 0.4, o\ = 90 k ms , and 02 = 500 kms 1 (Arzou- 
maman, Chernoff & Cordes[2002). For SBHs we adopt the same 
distribution, but with velocities smaller by a factor Mch/A/sBH 
(White & van Paradii JT99€l Gualandris et al l 20051) . The velocity 
kick leads to an increase of the eccentricity of binary and a change 
in its energy; the kicks are generally insufficient to ionise the binary 
systems (van den Heuvel et al l 200fJ) . 

The surviving binary loses orbital energy due to the emission 
of GWs. A CR of mass Mcr on an orbit with semi-major axis 
a and eccentricity e around an IBH of mass M, ^> Men loses 
energy due to the emission of GWs at rate of 
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If the orbit is circular (e = 0) spiral-in occurs on a time-scale of 
(Peter{T%3) 



morge ~ 256 G 3 Mi Mcr 
= 3.5 Gyr 
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For non-zero eccentricities the in-spiral time is shorter by a factor 

~(l-e) 7 / 2 . 

We perform binary population synthesis with the above de- 
scribed scenario to compute the fraction of high mass binaries that 
become potential LISA sources. The mass of the stellar compan- 
ion was selected from the initial mass function dN/dM* oc M^ a , 
where we assumed a = 2, consistent with the mass function in 
the core of clusters in which a runaway merger occurred (Portegies 
Zwart et al. I 20041) . We assume a minimal mass of IOMq, since 
lighter donors cannot account for UL X luminosities (Hopman et 
al. 2004; Portegies Zwart et al, 2004), and a maximal donor mass 
of 1OOM . 

LISA is likely to observe (IBH, CR) binaries in the phase 
before the merger, rather than the merger event itself. However, 
the CRs spent only a very short time in the LISA frequency band 
as compared to the Hubble time (see *|3j. The question whether 
these sources are detectable therefore depends on how many merge 
within a the Hubble time. 

In figures and J5J we show the results for an IBH of M, = 
3 x 10 3 M Q and M. = 1O 4 M , respectively. Only a fraction of the 
star spirals in within a Hubble time. For M. = 3 x 1O 3 M nearly 
all objects which spiral in are NSs, with only very few SBHs. For 
more massive IBHs, a significant fraction of SBHs also spirals in 
fast enough, and the total fraction of merging objects exceeds 10% 
for Mm ?L 10 4 M . In figure |5) we show /merge, the fraction of 
stars which spiral in within the age of the Universe, as a function 
of A/.. After the velocity kick, the CRs have eccentricities up to 
e < 0.9. While this decreases tmerge, these eccentricities are much 
smaller than the e ~ 0.995 eccentricities found by Hopman & 
Alexander (2005) for direct GW capture. By the time the stars spi- 
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Figure 1. Scatter diagram of the merger time t mcr gc as a function of the 
zero-age main-sequence mass of the donor A/zams- Events below the hor- 
izontal solid line merge within the age of the Universe. For this image we 
generated 5000 binaries with a a = 2 power-law initial mass function for 
the donor mass and we selected an IBH mass of M, = 3000 M@ . A small 
fraction (144 or 2.9%) of the objects generated experienced merger within 
a Hubble time. The sudden jump at 20 Mq indicates the transition between 
NS and SBH formation. 
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Figure 2. Same as figure (1), but with M. = 10 4 A/ o . In this case much 
more (~ 15%) stars spiral in within the age of the Universe, and in partic- 
ular there is a significant contribution from SBHs. 



ral in to orbital frequencies v > 10 4 s, which is in the LISA band, 



the orbits are close to circular. 



2.2 Interactions of the binary with cluster stars 

The two-body relaxation time t r of star clusters which host a run- 
away merger is short (t r ~ 30 Myr, Portegies Zwart & McMillan 
2002), and these evaporate on a time-scale of 10 s years if tidally 
limited. As long as the cluster has not yet evaporated the (IBH, CR) 
binary has interactions with cluster stars. Since the (IBH, CR) bi- 
nary is "hard" (Heeei d 1975ft . these interactions tend to decrease 
the orbital separation between the IBH and the CR (Mille d 20021) . 
In addition, scattering changes the eccentricity of the binary. The 
angular momentum vector J of the binary performs a random walk, 
and its magnitude J samples angular momenta < J < J m on the 



0.15 



0.05 



_j 

100 



I 000 



m bh [M S ] 

Figure 3. Fraction / me rge of CRs which spiral in due to GW emission 
within the age of the Universe as a function of the mass of the IBH. The 
error bars are estimated by the 1-<t Poissonian uncertainties of the simula- 
tions. 



cluster's relaxatio n time- scale (Alexander & Hopman 2003; Hop- 
man & Alexande rllobll) . Here J m is the maximum angular mo- 
mentum of a binary of given energy. When J decreases the orbit 
becomes more eccentric, and the in-spiral time decreases (see eq. 
0) 

It is not straightforward to quantify the effect of gravitational 
two-body scattering, since the cluster is not in a steady state, and in 
particular the relaxation time in the core can vary wildly. For sim- 
plicity we therefore neglect the effect of scattering in the follow- 
ing discussion on the in-spiral rate and the number of observable 
LISA sources. We note, however, that scattering may significantly 
increase the number of GW sources, so that the following LISA 
detection rates should be regarded as lower limits. 



3 OBSERVABLE GRAVITATIONAL WAVES FROM IBHS 

The dimensionless strain of the GWs emitted at a frequency v — 



nr 
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s 1 from source at a distance d — dm-pc Mpc is 
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(e.g. Sigurdsson & Ree d 19971) . LISA is sensitive to frequencies in 
the range 10" 4 Hz < v < 1 Hz. At v = 10~ 3 Hz, LISA can detect 
sources with strains larger than h = ft-23l0~ 23 , where /i—23 ~ 1. 
This estimate is based on a 1 yr observation with signal-to-noise 
ratio S/N = 1 (see, e.g., http://www.srl.caltech.edu/lisai. In the fol- 
lowing we assume for concreteness a GW source with frequency 
v = 10~ 3 Hz; application to other frequencies is straightforward. 
Sources can be observed to distances up to 
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We assume that only ULXs with lumi nositie s Lx > 
10 40 ergs -1 contain IBHs (Portegies Zwart et al l 20041) . Presently 
the average number of > 10 40 ergs -1 ULXs per galaxy is ~ 0.1 
(Swartz et al j 20041) . The star formation rate dropped by an order 
of magnitude since z ~ 2 (Madau, Pozzetti & Dickinson l 19981) . 
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ULXs correlate with star formation; here we assume that the num- 
ber of ULXs is proportional to the star formation rate, in which case 
the number of luminous ULXs per galaxy at an earlier time was 
Njjlx ~ 1. Since the ULX lives for a time t* = 10 7 yrt7, this 
yields an ULX formation rate per galaxy of 1CP 7 yr~ 1 A r uLxif 1 - 
As was discussed in the previous section, only a small fraction 
/merge = 0.1/_ i of the ULXs with M. < 1O 4 M leave behind a 
remnant binary which spirals in within a Hubble time. The rate at 
which observable GW sources are produced per relic ULX is then 

r GW = 10 -8 yr-7-iiVuLxtf 1 . (8) 

At the point where the period P = 2na 3/2 /(GM.) 1/2 equals 
10 3 s, the in-spiral time is 

M. 



th = 155 yr 
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where we assumed that the orbit has circularised when the fre- 
quency is this high. We thus find that the mean number of sources 
emitting GWs with frequencies of v ~ 10 _3 s _1 per galaxy is 
M\ = TawtL, or 
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The local galaxy densit y is es timated to be n ga i w 3 x 
10 -2 Mpc -3 (Marinoni et al.Q999). If the maximal distance at 
which the GW source can be observed is d ma x, the number of 
LISA sources in the sky at any moment is given by A/tot = 
jVin ga i47rd max /3, or, using equation Q, 



harden the binary, and change its eccentricity. As a result t me rge de- 
creases, and thus the number of potential LISA sources increases. It 
is not implausible that nearly all (IBH, CR) binaries merge within 
a Hubble time, in which case /_i = 10. Second, the life-time of 
ULXs is probably significantly shorter than the main sequence life- 
time t* of the star. A more realistic assumption would be that ULX 
luminosities are only achieved when the donor is near the terminal 
age main sequence, in which case ty ~ 0.1, boosting the predicted 
detection rate for LISA with an order of magnitude. In conclusion, 
the number of observable GW sources could easily be orders of 
magnitude larger than the expression WW indicates. 

Previous studies of the number of potential LISA sources from 
binaries with an extremely small mass ratio focused mainly on 
cases in which the orbital energy is dissipated by t he GW s them- 
selves (e.g. H ils & B endeJ 1 99 a Si gurdss on & Ree d 19971 Ivanov 

E 2002 1 : Miller ! 2002t Freitag l 200 llj 20031: Alexander & Hopman 
2003; Hopman & Alexander 2005). In that case the event rate is 
on the order of a few per Gyr, and the orbits are highly eccentric, 
with typical eccentricities as large as e ~ 0.995 for v = 10~ 4 s -1 
near an IBH of M. = 10 3 M Q (Hopman & Alexande j2005l) . In 
that case stars spiral in very quickly, and emit GWs in the LISA 
band only for a short time, in contrast to the here discussed tidal 
capture scenario. 

GW sources originating from tidal capture sources have nearly 
circular orbits when they enter the LISA frequency band, leading to 
a distinctly different signal than the highly ecce ntric s ources origi- 
nating from direct capture (e.g. Barack & Cutlei2004; Wen & Gair 
2005). LISA can determine both the mass of the IBH and the ec- 
centricity of the orbit. Detection of GWs from an IBH will give 
proof of the existence of these objects. If the signal stems from an 
orbit with low eccentricity, this supports the scenario that ULXs are 
accreting IBHs in binary systems. 



4 DISCUSSION 

For IBHs of M. < 3 x 1O 3 M , /_j ~ 0.3 and almost no SBHs 
spiral in (figure (3). The main contribution to LISA sources comes 
from NSs, i.e. Mcr/IOMq ~ 0.14.We estimate that the merger 
rate for NSs with IBHs is too low to be likely to be seen by LISA. 
However, in this estimate we ignored the effect of three-body scat- 
tering on the binary orbit. Though hard to quantise, its result may 
be a marginal detection rate for NSs. We also note that we esti- 
mated that SBHs have a velocity kick distribution similar to that of 
NSs, but with velocities smaller by the mass ratio Mch/^fsBH of 
the objects. This estimate is rather uncertain, and may be too con- 
servative. Larger kick velocities cause more SBHs to merge within 
a Hubble time, in which case the contribution of M, < 3 x 1O 3 M 
IBHs to LISA would increase. 

For IBHs of M. > 3 x 1O 3 M we find that /_i ~ 1 and 
for SBHs tmerge decreases to well within a Hubble time. In this 
case stellar mass black holes (Mcr/1OM ~ 1) give the most 
promising GW sources, with a detection rate of about 10 per year, 
in the case that three-body scattering is ignored. Including three- 
body scattering in this estimate may boost the detection rate with a 
sizable fraction. 

Our estimate in equation for the number of sources that 
LISA can detect is conservative. First, in Sj3]the (IBH, CR) binary is 
considered to be isolated. Indeed the host cluster eventually evap- 
orates, but this is preceded by a phase during which the binary in- 
teracts with other cluster stars ( i!2.2> . These interactions tend to 
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